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Abstract 


at ONERA, which makes up the p%load df fhe D5B Caltor sS 

French Space Agency) confirmed the existence of an efectric curj^f char4^^^^ 

proof mass under the influence of the magnetospheric protons and revenlfd f 

from dat. oB proton (lin. at the aatelllte altSsTxhe coSSarSor.^^ 

fo f»hows that the calculation method is valid and precise enoueh 

to be used for drag-free or accelerometric satellites. precise enough 


I. INTROIrt’CtldN 

A "drag-free" satellite^ is essentially made of a probf mass protected from 
the surface forces acting on It and which are due to the slowing down created by 
the residual atmosphere as well as the various radiation pressures. The piloting 
system of the satellite controls the thrusters in such a way that the cage containing 
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the proof mail never comoi Into contact with it. The iatolllte trajectory is then 
In principle purely gravitational. In practice, the componiutlon by the thruitei i 
of the iurface forcoi ii not perfect, the reilduc oi compeniatlon being aue to 
internal perturbatlng forcei acting on the proof matis and thus modifying the 
reference trajectory. 

The itrueture of an accelerometrlc .satellite is .similar, but the servo-eontro 
function takes place by the action of a force on the proof ma.s.s so that the latter is 
maintained in the vicinity of the cage center. The measurement of the internal 
forces developed by the servo-dontrol then constitutes a measure of the sum of the 
surface forces acting on the satellite. The internal perturbatlng forces are also 
at the origin of the physical limitations pertaining to this type of instrument. 

Electrification of the proof mass by proton and electron fluxes from the radia- 
tion belts makes up one of the main perturbations, all the more so as this electrifi- 
cation may increase with time and reach high levels. 

So, during the definition of a drag-free or accelerometrlc satellite, <t is 
important to be able to determine a priori the value that will take the current 
charging the proof mass in orbit so as to decide on the procedure to implement to 
compensate this effect. 

The present paper gives the results of a comparison which has been made 
between this charging current as calculated for a particular orbital configuration, 
and the results of measurements obtained in orbit on a three-axis accelerometer 
(Cactus) making up the payload of the French satellite D5B-Castor. 

2. OF TIIF nitXii-FHFF \M) XCt FI.FIlOMFl KH, SXTFI.I.H F.< 

Let us consider (Figure 1) a material sphere of mass m placed inside a cage 
fixed within a satellite. The mass of this satellite - including that of the proof 
macs - is M. The position of the proof mass center Og is defined by the vector 
I in a reference frame OgXYZ linked io the satellite and such that Og be at the 
center of mass of the satellite. 

Let It be: 

F the resultant of the internal forces of attraction of the proof mass by the 

La 

satellite, 

Fg the resultant of the surface forces acting on the satelllt- (atmospheric 
drag, radlatlou pressure), 

Fp the thrust due to the thrusters, and 

Gt, and G , the local gravitational acceleration in Og and Og. 



If and represent respectively the coordinates of Op and Og in an ab- 
solute frame of reference, the movements of the proof mass (m) and the satellite 
alone (Ri - m) are given by the equation of dynainics: 




" M - m 




M 


( 1 ) 


( 2 ) 


2.1 Ihrap'i'rec SulrIliU* 

If Vve suppress any link between proof mass and cage (F^ - 0), the trajecto.'.'y 
of the proof mass is purely gravitational. 

By piloting the satellite in such a way that, under the action of the thrusters, 
the amplitude of T remains at any moment lower than a value previously 
chosen, we have a satellite whose trajectory can also be characterized as purely 
gravitational, as it only differs from that of the proof mass by a distance almost 
equal to a distance always small as compared to the satellite dimensions. 

But if a perturbatlng force T remains between the two bodies, the true trajec- 
tory of the satellite departs from a purely gravitational trajectory and the metric 
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■ 2 Xi-rrli-rnnh'IrM 

u< ‘i '• '.'lU’lHtp is of any propulsi<.m moans (1''^^ 0), but a Horvo-ronti'ol 

iH. ,Mu.n Mcts nil the proof mass to maintain the points and in a sin^^le 

|. w.Uinn (\*iM in the erro.^ of the servo-control). In general, considering the small- 
n* > - ..f tho int i'os, no material contact should exist between proof mass and cage 
;i i ti ♦* liaisoTi force F is — for example —of electrostatic nature. 

In »him- mnditions ^ 0 and Gj^ G^. 

linlntir-Tv^ (1) and (2) then give: 


M 


(4) 


j'lin 'U (‘cle?'ntinn imposed on the satellite by the external forces — apart from 
i. itv - i- c(pial to the force of proof mass-cage liaison divided by the mass of tho 
•< * rnrts-. 'j'l.iis, the measurement of this liaison force Fj makes it p(;ssible to 


• ’« Mio resultant of external forces Fj,, 

1*' :>n intci nal pei'turbating foi’cc f is added to the force Fj developed by the 
r\. -« ontrol^ the measure of the acceleration due to the external forces is then 
‘1* iil-e by a systematic error equal to f/m. 


t'i HI I Him ION III i Ml nil: i i.i:( itui frxnoN ot iiik iMiooi 

It) I'oMj systems that have been just described, the perturbation (of the trajec- 
• r \ '*r f>r the for* e measurement) is directly given by the acceleration f/m that 
li* pert nrl a ting force would communicate to this proof mass alone. 

ilie-f pot turl)riting forces are of various natures and have already been the 

3 * 5 fi 

I i . t ,.r fF tailed studies as well as of measurements in flight.*** * Among them 

^ • . rr provided by the defect of electric neutrality of the proof mass constitutes 

f *!iirbation that may become very important. In this f'ase, the proof mass is 

i |ortfd to an attractihg force by the cage walls on which are induced electric 

barges wliose sum is equal and of contrary sign to the charge carried by the proof 
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maa.s. Due to the .spherical .symmetry of the eaf?e, thi.s attraction i.' a central 
force around a point O, (Figure 1), geometric center of the cage which we try to 
make a.s neat as po.s.stble of center 0^.. This periurhating force f,, can he ex- 
pre^iHed byj^' ^ 

fjj mllQ^tr -X) 

where T t)^)j and ji is a coefficient defined hy the geometry of the instrument. 

This force is thus proportional to the square lif tiie eiectrie charge carried b> 
the proof mass. 

This may have two very dirferoftt origins. 

(1) Electrification of internal origin which appears when the proof mass leaves 
its contact with t.he cage, a contact obtained eiti»er in the presence of gravity (on 
the ground) or under the action of a force obtained by remote control (in orbit). 

This electrification is due to the charges developed either by instantaneous poten- 
tials of the various electrodes of the cage or by the differences of the work func- 
tions of the materials making up the proof mass and the cage walls. 

(2) Electrlficaticvi of external origin due to the accumulation of charges pene- 
trating into the satellite, originated by the high energy particles of the radiation 
bells. 

The first kiftd of electrification may be minimized by appropriate toclmoiogioal 
means. The second appears as a current charging the proof mass witose order of 
magnitude is hardly predictable without a detailed study. Indeed, it i.s very ditfi- 
cult to know a priori if the proof mass charge will reach a prohibitive value within 
a few* days or few years. 

Though it is necessary, during the project of a drag-free satellite or an 
accelerometric satellite, to foresee the adequate means for discharging the proof 
mass, it is highly desirable that their optimization might take into accoimt the 
maximum and minimum values of the electric current that will charge the proof 
mass III orbit. 


I. TiiKCu ii .' n.Mi 


The Cactus accelerometer (in French; Capteur accelerom^trique capacitif 
triaxial ultra sensible) has been designed and built by ONERA, and made up the 
payload of the French satellite Castor (D5H) placed in orbit on 17 May l!'7ri. T1 
satellite has been built and launched by CNES (The French .Space Ageney) who 
ensured the further ekploltatlon of the instrumentation. 
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1,1 — Ui‘*>rri|ilMm ol’ ilti* 

The Cactus accelOhometer has a measuring range of tlO * g on each oi its 
three axes. The sum of intcinal perturbations has been eVi9luated before latnch- 

•9 

ing at 10 g. These values, as well as all the other characteristics, have been 

7 

confirmed by the results obtained in orbit. * • 

The core of the accelerometer is made by a proof mass in rhodiated platinuin 
whose mass is 550 g — placed in a cage forming with it a gap of fi5 ^m (Figure 2). 
The force linking proof mass in cage is of electrostatic nature and is obtained by 
means of continuous voltages applied on three systems of electrodes distributed 
over three orthogonal axis. These voltages are made proportional to Ih^ relative 
displacement of the proof mass in the cage thanks to another set of elector des 
realizing, on each axis, a capacitive nu?asurement of position. The system func- 
tions by position servo-coritrol of the proof mass and the measure of the voltages 
applied on the acting electrodes on each axis makes it possible, after preliminary 
calibration, to know the liaison force and thus to determine the sum of the 
external forces (Eq. (4)). 



f* l;^re 2. Cage, Proof Mass and Electrodes of 
the Cactus Accelerometer 

The D5B/Cactus experiment aimed at: 

(1) ensuring qualification in orbital flight of the accelerometer, ^ and 

(2) providing scientific data on aerunomy. 

The orbit chosen was slightly excentric, with an inclination of The alti- 

tudes of apogee and perigee of the first orbit were respectively 1275 km and 277 km 
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1,2 (iiilralaiiunH of tho I'roof Mohn TUrronl 

At the same time as the studies for defining this experiment, theoretical and 
experimental studies^* wefe performed with a view to attempt to determine 

the mean value of the proof mass charging current. 

These Wot*ks showed that; 

(1) the evolution of the proof mass charge is essentially due to the bombard- 
ment of the satellite by the magnetosphere electrons and protons; and 

(2) the interactions of these two types of particles with the satellite structure 
have w^ldely diffet*eht characteristics. 

4.2.1 ELECTRONS 

While crossing the materials of the satellite and the accelerometer cage, the 
primary electrons give rise to secondary electrons and to photons,, part of which 
reach the proof maSS artd interact with it. The results are that electrons and 
photons circulate in both directions between proof mass and cage. The opposite 
fluxes of electrons are not equal, hence the existence if a charging current whose 
value and sigh can be determined only by a detailed study. 

4.2.2 PROTONS 

Contrary to electrons, the protons crossing the matter do not generate second- 
ary effects of any importance, and propagate practically in straight lines. The 
protons stopped within the proof mass are at the origin of a charge increase. 

The work carried out at ONERA showed that; 

(1) the charging current is essentially due to the primary effect of the protons 

stopped by the proof mass; the presumed mean value of this current has been eval- 
uated at +(2 ± 1) Coulomb per day; and 

(2) only the electrons whose energy is about 4 MeV can give a perceptable 
charging current; by extrapolation above 4 MeV of the known values of the flux, 

the presumed electronic charging current has been evaluated at -1,7 10 CoUlomb 
per day, as a mean value. 

At the end of this study, the expected mean value of the charging current was 
thus near f2. 10“^^ Coulomb per day, or 2.3. 10*^^ ampere. 

The methods used to calculate this current, ns well as the assumption adopted, 
are described in Section 5. 

i.3 ilHermimiiiurt in Urbilal of the Pro<if Muss rhur^sinp (Current 

The acceleration measured by the accelerometer when the satellite is near its 
apogee — where the atmospheric drag Is negligible — and In the shadow of the earth — 
where the acceleration due to the sun radiation pressure disappears — constitutes a 
good measure of all the Internal perturbations of the instrument, as the only error 
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of this measurement comes from the earth radlatloh pressure, which provides the 

satellite with an acceleration of the order of 4. 10 g. 

the existence of a proof maSs Charging current is well revealed by an increase 
with time of this acceleration and by the returns of the latter to its bottom level 
during each contact between proof mass and Cage, obtained by remote control. 

But, moreover, systematic readings of these data have also revealed time periods 
of about . 10 days, renewed every 38 days, and during which the charging current 
betdm&s weaker. 

As an illustration. Figure 3 represents the values of the modulus of the accel- 
eration measured by the accelerometer when the satellite is at a high altitude 
between the 10th of M&y and the 20th of August 19t6. We can see on this figure the 
periods when this current weakens: they are the periods from 18th to 31st of May, 
from 25th of June to 8th of July and from 28th of July to lOth of August. Outside 
these, the charging current takes again a higher value characterised by the increase 
of the measured acceleration. 



Figure 3. Acceleration Measured by Cactus at High Altitude 


Various methods have been used to determine the charging current from the 
data transmitted by the satellite.'’* “ These methods consist in identifying the 
various Internal perturbations and the accelerations due to the external forces 
from a realistic modeling of these accelerations and by using the attitude data of 
the satellite. When the electric charge level becomes high enough, we consider 
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that th6 pertOhUtlort dufi to the charge constitutes the hjain term of the accHcrntion 
measured when the satellite is at a high altitude. I'sing Kq. ('•,) on the smoothed 
data theft permits a simple but sufficiently precise calculation of Q and its time 
variation, for the appUcation of Eq. (5), the knowledge of tiie values of /{ and T 
IS necessary. Coefficient H had been determined on the ground beP.re launch and 
has a value of B.9. 10 A s . The compOhents of the T vector have been meas- 
ured in flight by a particular manoeuvre'’; its modulus has the value 2. Hi am. 

These various methods gave coherent results which are as follows; expressed 
as a mean value of the charging curj-ent during the considered period of ti.ne: 

• For periods with strong char *g iiig current : 

. From 28th of June to 8th of July 1975: 2. 1. lt)-“ Coulomb per day (2.4. 10-*'’A) 
.From 3rd to f.th of November 1975; 1.3.10’“ Coulomb per day (1.5. A) 

.From lOth to 20th of June 197r 1.3.10’“ Coulomb per day (l‘. 5. 10’^^’ A) 
.From 25 August to 5 September 1976: 1.02. 10’“ Coulomb per day 
(1.2. lO’^'^ A). 


- For periods w ith weak charging current: 

.From 23 rd to 28th of June 1976: current lower than 5. l0’“ Coulomb per dav 
(5.8. 10’^“ A) 

. From 6th to 1 5th of September 1976: 1. 2. 10’“’ Coulomb per dav 
(1.4. 10 ’“a). 


.*>. iiKr\ii.i:i) i vi.ci i.vnoN nii; imuhh m\s> t ii\n(,i\(, u mu;M 


3.1 lillerpn-tuliun of (ho Xuriulion^. Oll^.l■rvl>d on th.- (ihar^iii^ Cum-nt 

We can see that the periods during which the charging current is strongly 
attenuated are centered on the dates when the apogee latitude i North and at its 
maximum value, that is 30°. a value corresponding to the orbit inclination. This 
important variation of the charging current may be explained by the following fact; 

(1) the proton and electron fluxes decreasing rapidly with altitude, the charg- 
ing current reaches a significant value only when the satellite Is around its apogee; 
and 

(2) due to the fact that the magnetic anomaly of South Atlantic which is char- 
acterized, at the altitudes where flies the satellite, by more intense panicle fluxes 
centered oirer a point situated at about 25° latitude south and 40° longitude west 

(Figure 4). the charging cuhrent takes a high value each time t.he satellite Hight 
crosses this zoae. 

Thus we can see that when the apogee latitude is around 30° south the satellite 
crosses this zone every day, while when the apogee is about 30° north this zone is 
avoided by the satellite. 


I 


I 



Figure 4. Proton Flux Contours -E > 100 MeV 


This effect, which is a combination of the orbital movemfent with a^ gcographic 
inonialy, .s all the more marked as the orbit inclination is close to 25^ which is 
thr* rase of the Castor satellite. 

With a view to provide a quantitative support to this interpretation, studies 
in ovinusly carried out on the determination of the charging current have been 
. vpcoted in order to evaluate the daily charge acquired by the proof mass during 
-omplete -cycle cf 38 days. 

i.J i ult iiliilioH of \hv h> dir >alrllilr 

Since 1072. when the first calculations of the charging current have been made, 
the description of the radiation belt has been refined and the extrapolation of the 
calues of electron Huxes at energies higher than 4 MeV (see Section 4. 2) appeared 
n. having no meaning. This statement confirms that only the protons in period of 
•i..t !o:il activity are at the origin of the charging current. 

Two methods a *e usually used for evaluating a flux received by a satellite 
during its useful life. If the mission duration is long enough, the experi- 

cuter may he interested only in the mean flux and me calculation consists in 
determining the probability for the satellite to pass within each volume element, 

I, delated m energy. Here the calculation requires more precision. Indeed, a 
satellite in low orbit, as D5B. i^ subjected to an inten„e particle bombardment 
,.nlv during rather short periods which correspond to the passage through the Sov.th 

AHrintir nnov>aly. 

1 u cab ulnte the proof mass charging current, we must bo able in each point 
. t -he ..t bit to ostimate the Hux of incident protons. Account being taken of tlie 
l, ,o ..f mbital parameters, asccndln- node and perigee argument, a calculation of 
,ho number of protons received every 24 hr seems sufficient and •.emains signif- 
„ .mt To this end, each orbit is described step by step and every minute tho 
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fii‘ographic roordlnatcH (altitude, latitude, InnMitudel are transfoi nu d into 
lua^rnetic coordinates <Il, Id that permit the consultation of tije files of protrin-> 
whose energiC».s are higher lOO MeV (see Section fi, :0. 

Figure 5 gives as a function of time, frc>m :trd of June to :trd of Julv, tlu* 
omnidirectional mean fluxes of protons of an energy higher than lOtt MeV, per 
square centimeter and peh day, as well as the apogee latitude. The first curve 
illustrates well the rtorth-south asymmetry of the radiation belt due tn the South 
Atlantic anomaly at the altitude considere d. The poritKlicity of this phenomenon 
makes it possible to extrapolate this flux curve before .'ii d rtf June and after ii d 
of July, as represented by a broken line the figure. 



Figure 5, Isotropic Proton Flux Density Heceived by 
the Castor Satellite and Apogee I.atitudc at the Same 
Periods 


o..‘i UiM»r|i(inti <if l*rotoh> liy tlir >ulelliti* Mrm turr 

The satellite shape is a regular polyhedron with 2t» faces whose geometric 
center is at the center of mass and also at the center of the proof mass. A simple 
model of proton absorption by the satellite structure has thus been established by 
decomposing it into 2(i equal solid angles. For each of thcsi\ the nature and the 
thickness of the vaHous materials encountered by a particle moving on a straight 
line and reaching the center have been surveyed. Fsing the density of each matei ial 
as a weighting parameter, these various thicknesses have then been converted in 
equivalent thicknesses of aluminum. These equivalent thicknesses vary fmm 270 
mm (in a solid angle of 4:r/l.'i) to 42 mm (in a solid angle of 4 t/20), 
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Moreover, the accelerometer proof mar»s beih^» in platinum and liavin^ a diam 
eter of 3U mm, thO eq»d»^alent aluminum thickncris tUirre.^pondinf; to a diametral 
cru^sin^' of the proof nia-s^ is .'1U» mm. 

Consequently, the proof mass char^'t? will he due tc» proton^ ener/»y is 

hi^^h enough to cross 42 mm of aluminum but whose ener r^y remains lowei’ than 
that necessary to cros i rjlbl mm of the same metal. 

Figure (», taken from, makes it prjssible to give a simple analytic formula- 
tion of the path x of protons of energy F in aluminum; 

X o E 

With X in mm and E in Me% we obtain the following empirical values: 
o 0. 01 and ^ 1. 73. 

Expression (H) thus permits for each solid angle of the satellite the determina 
tion of energy band of the protons which will participate in the proof mass charge. 
For the satellite as a whole, there are the protons whose energy is between 100 
and 500 IVleV which are to be Considered. 



Figure t». Proton Path in Aluminum 


i t!ah ulalion »f tin* Ihiilv nmrjiin^r «f ihr l^fiad 

The argument is that used by J. Tiffon, recalled hereafter. 

The flux density of protons of an energy higher than E, o< »E), is ! ep»’(‘sente<l 
by the following analytic expression: 


IRC 


(7) 
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Parainotori k and i] arr doiorniinrd <‘m|iirir‘aUy from fhr hn<* wkieh 

nmko;s up a appivixiinaliun of iho lavs ivpixvsunfin^* tlio variation of Uu* lo^*a- 
rithni of tho inti*^r*att‘d Uux a s ' iunoti >>n of tJu* lo/[;aritiwii of <»nc*r^y. 

Ah an illuritr*ation, l i^^irc* 7 roju*osonts flu* valuer; of fho intf*^r'atod Ilux of 
protoftH rtHH ivod by Uio satollito durtiH* iho riay of atli of June liiTii. a w( J1 as the 
■straight line (.dvin^f its apprfjximate i XfU'Cs ;ion. 

Inside a solid aiud^^ *..*j (ld^»ure «) t oi ta*sf)nndin^» lu thi* previtjusly defined 
dividing? and h»r whicli the aluminum thiekness to bv ( vuhmh] is \ , the number t»f 
protons issued from a solid an^le dw and i eaehin^ the proof mass within a day is 
defined from Kq, ((J) and (7) b\ ; 






(8) 


whfere H is the mean radius of the screen of thickness Xj. AnumK the.-,e. a certain 
quantity dSi.,^ eomes out of the proof mass: 


dNi;, k 
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Hl’j cos t) . 
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’n this expression, r is the proof mass radius - in platinum - and x its 
equivalent aluminum thickness. ^ 

A number of protons stopped daily by the proof mass and penetrating into the 


solid angle s.v is thus: 
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or, after calculation: 
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Figure 7. Proton Flux Received on 
June 8, I07r; 


^ A 



Figure 8. Method for Calculating the Proof 
Mass Charge 


The application of Eq. (10) to each of the solid angles defined in Section 5.3 
and for each day during which the radiation doi»es received by the satellite arc 
characterized by the parameters k and n. permits, after summation, the calcula- 
tion of the electric charge acquired dally by the proof mass. 

These results are presented on Figure P, which represents the daily charge 
acquired as a function of time. The abscissa corresponds to the 30 days of the 
month of June 1070 for which the calculation has been performed. The dates 
where the apogee latitude is at 30^^ north and 30^ south have also been shown. This 
makes it possible to generalize the result obtained at any date of the satellite life 
in as much as the orbit decay is hot too important, 
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The valnee of the charging eurrente deterntined fror, the data provided hv th 
aeee lero„eler are aleo ahown on Figure 0 aa hortrontal llnea. Thlae ilnea glvl 
mean value of thla current over a period r„rrea,va„dlng to the line length Tl 
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rurient as a funrtion of orbttal situations. ^ 





i 


Thin evaluation should permit a better definition of the means to imple.oen. 
and the procedure to use to malrttaln the proof mass electrification to a tolerat-b 

level. 
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